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. We apply the non-equilibrium fluctuation theorems developed in the statistical 

physics to the thermodynamics of black hole horizons. In particular, we consider 
a scalar field in a black hole background. The system of the scalar field behaves 
. stochastically due to the absorption of energy into the black hole and emission 

5h ■ of the Hawking radiation from the black hole horizon. We derive the stochastic 

" " " equations, i.e. Langevin and Fokker-Planck equations for a scalar field in a black 

hole background in the /i — )• limit with the Hawking temperature UkI^tx fixed. 
We consider two cases, one confined in a box with a black hole at the center and 
the other in contact with a heat bath with temperature different from the Hawking 
temperature. In the first case, the system eventually becomes equilibrium with the 
Hawking temperature while in the second case there is an energy flow between the 
black hole and the heat bath. Applying the fluctuation theorems to these cases, 
we derive the generalized second law of black hole thermodynamics. In the present 
paper, we treat the black hole as a constant background geometry. 
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Since the paper is also aimed to connect two different areas of physics, non- 
equilibrium physics and black holes physics, we include pedagogical reviews on the 
stochastic approaches to the non-equilibrium fluctuation theorems and some basics 
of black holes physics. 
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1 Introduction 



The analogy of the space-time with horizons and thermodynamic systems have been 
extensively investigated, especially, in the black hole thermodynamics [1]. A black hole 
behaves like a blackbody with the Hawking temperature Th — hK/2n [2], and energy 
flowing into the black hole can be identified as the entropy increase of the black hole. 
Here, k is the surface gravity at the horizon and the entropy of the black hole Sbh is 
proportional to the area of the horizon A as Sbh = A/4G in the Einstein-Hilbert theory 
of gravity. The thermal behavior is essentially quantum mechanical. Furthermore, such 
a thermal behavior is not restricted to globally-defined horizons like an event horizon 
of a black hole, but also applicable to local horizons such as the Rindler horizon of a 
uniformly accelerated observer. At the quantum level, the notion of horizon entropy 
must have more fundamental meanings, since it gives transition rates of area-changing 
irreversible processes of black holes, and will be related to the quantum statistical nature 
of the space-time. Such microscopic views have been proposed in string theory, especially 
in the approach based on D-brane constructions [3]. 

Because of the thermodynamic behavior of black holes, a system in a black hole back- 
ground behaves as a system in contact with a thermal bath. In particular, if we consider 
a (scalar) field in a black hole background, its effective equation must be described by a 
stochastic equation with dissipation and quantum noise. The dissipation comes from the 
classical causal property of the horizon; the black hole horizon absorbs matter and, once 
they fall in, they cannot come out. The property is the basis of the membrane paradigm 
of the black hole [4], in which Ohm's law or the Navior Stokes equations hold on the 
membrane at the (stretched) horizon. On the other hand, the noise term (or fiuctuation) 
comes from the Hawking radiation, which is essentially quantum mechanical and, hence, 
we need to quantize the system in the black hole background in an appropriate way. The 
first purpose of the present paper is to derive such a stochastic equation of motion for a 
scalar field in a black hole background. The stochastic equation of motion of a string is 
previously derived in [5, 6] based on physical intuition of the Hawking radiation, or in [7] 
by using an analogy with the Schwinger-Keldysh formahsm in the context of AdS/CFT 
correspondence [8]. Our approach is similar to them, but we obtain the effective equation 
by explicitly integrating fluctuating degrees of freedom. Namely, we introduce infinitely 
many variables between the horizon and the stretched horizon and consider them as envi- 
ronmental variables. By integrating them, we can show that the variable at the stretched 
horizon behaves stochastically with a noise term. Though the environmental variables 
are living outside of the horizon, they can encode information in the black hole through 
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choosing the Kruskal vacuum with the regularity condition at the horizon. In this sense, 
the integration of the environmental variables corresponds to integrating hidden variables 
in the horizon. The derivation of the Langevin equation is one of our main results. 

The second purpose of the paper is to apply the non-equilibrium fluctuation theorem 
[9]- [10] developed in the statistical physics to the scalar field in the black hole background. 
In thermodynamic systems, entropy is always increasing (or remaining a constant). But 
for a mesoscopic system where fluctuations are large, there are nonzero probabilities that 
the entropy of the system decreases. The fluctuation theorem relates probabilities of en- 
tropy decreasing processes to those of entropy increasing ones in terms of the equilibrium 
thermodynamic quantities. It is a very general theorem that can hold for various dynam- 
ical and non-equilibrium systems including classical Hamilton dynamics in contact with 
a heat bath, stochastic equations with dissipation and noise, or quantum mechanical sys- 
tems. The Jarzynski equality can be derived from the fluctuation theorem, and the second 
law of thermodynamics is implied from the Jarzynski equality. We use the word implied 
here because the second law can be derived only if we assume that a system is relaxed 
to an equilibrium state after a long time. An application of the fluctuation theorem to a 
scalar field in a black hole background is straightforward once we obtain a stochastic equa- 
tion of motion. We can derive the generalized second law of black hole thermodynamics, 
or Green-Kubo formula of the linear response and its nonlinear generalizations. 

The paper is organized as follows. In section 2, we briefly review the stochastic ap- 
proach to thermodynamic systems, Langevin equation and Fokker-Planck equation. An 
important property of the stochastic equation is that it violates the time reversal sym- 
metry which can be measured by an entropy increase in the path integral. In the next 
section 3, the fluctuation theorem for a stochastic system is reviewed. It relates the en- 
tropy increasing and decreasing probabilities. From the fluctuation theorem, the Jarzynski 
equality is derived. In section 4, we derive an effective stochastic equation of a scalar fleld 
in a black hole background. In deriving the Langevin equation, the quantum property of 
the vacuum with a regularity condition at the horizon is very important, which is first 
explained. We then introduce a set of discretized equations of a scalar field near the black 
hole horizon, and integrate the variables between the horizon and the stretched horizon. 
The integration leads to an effective stochastic equation for a variable at the stretched 
horizon. This has the same spirit as deriving a Langevin equation of a system in contact 
with a thermal bath [11, 12, 13]. In section 5, we apply the fluctuation theorem to the 
scalar field in a black hole background. We consider two different situations. In the first 
case, we put the scalar field and the black hole in a box with an insulating wall. By 
applying the fiuctuation theorem, we can derive a relation connecting entropy decreasing 
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probabilities with increasing ones. The ratio is given by the difference of free energies. 
Prom this, the generahzed second law of black hole thermodynamics can be derived. In 
the second case, the wall is assumed to be in contact with a thermal bath of a differ- 
ent temperature which is slightly lower than the Hawking temperature of the black hole. 
Then there is an energy flow from the black hole to the wall. By applying the fluctuation 
theorem to it, a linear response theorem of an energy flow to the temperature difference 
can be obtained. In the appendix A, we review a derivation of the path integral form of 
the Fokker-Planck equation. In the appendix B, we will discuss the relation between the 
noise correlation and the flux of the Hawking radiation. In the appendix C, we explain 
the fluctuation theorem for a steady state and derivations of nonlinear generalizations of 
Green-Kubo formula. 



2 Stochastic Equations of Motion 

We first briefly review stochastic approaches to classical statistical systems. In particular, 
we focus on the path-integral representation (Onsager-Machlup formahsm) of the Fokker- 
Planck equation and emphasize the role of time-reversal symmetry. Readers familiar with 
non-equilibrium statistical physics can skip this and the next sections. 



2.1 The Langevin Equation 

The Langevin equation is a phenomenological equation of motion of a particle with a 
friction term and thermal noise. It is commonly described as 

dV 

mv^-^v- — + ^. (2.1) 

V{x) is an external potential for the particle. 7 is the friction coefficient and ^(t) is 
a thermal noise (or a random force) which is often assumed to have a Gaussian and 
white-noise (delta-correlated) distribution 

(C(t)) = 0, mi{t')) = '^lT5(t-t'). (2.2) 

The coefficient is determined to satisfy the equipartition theorem with the temper- 
ature T through the fluctuation-dissipation theorem. The noise average (•••) can be 
represented by the following path integral 



{Fit)) = J V^F{t)exp [-1 J dUdU{ti)^-^-^at2) 



(2.3) 
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with a normalization condition (1) = 1. If necessary, we can easily generalize the noise 
correlation to an arbitrary colored non-Gaussian noise. An well-known example that can 
be conveniently described by the Langevin equation is the Brownian motion of a particle 
or thermal fluctuations of an electric circuit voltage. 

2.2 The Fokker-Planck Equation 

From the Langevin equation, we can derive another type of a stochastic equation, the 
Fokker-Planck equation. It describes a dynamical evolution of the probability distribution 
P{X,t) of observables X at time t. Here X represents the variables {x,v = x). If the 
process is Markovian, i.e. the next state is determined only by the present state, the time 
evolution of P is given by the following Master equation, 

dtP{X, t\Xo, 0) = y" dX' [w{X' X)P{X', t\Xo, 0) - w{X X')P{X, t\Xo, 0)] . (2.4) 

Here P{X, t\Xo,0) is a conditional probabihty to find an event X{t) — X that has started 
from the initial value X{0) = at t = 0, i.e. P{X, t = 0\Xo, 0) = 5{X-Xo). w{X' X) 
is a transition rate from X' to X, which can be related to the Langevin equation in the 
following way. The first and the second terms of the right hand side of eq.(2.4) describe 
an incoming and outgoing fluxes of X respectively. 

The Master equation can be brought into the Kramers- Moyal form as 

dtP{X,t\Xo,0) 

^- j dr [w{X ^ X + r)P(X, t|Xo, 0) - w{X - r ^ X)P{X - r, t\Xo, 0)] 
^- I dr[l- e""^^] w{X ^X + r)P{X, t\Xo, 0) 

= E -^9x [CniX)P{X, t\Xo, 0)] , (2.5) 

n=l 

where we have defined 

Cn{X) = J drr^w{X^X + r)^ hm^^^{{X{t + At)-X{t)r)\xit)=x. (2.6) 

In the last line, we have rewritten the n-th moment of the transition rate by a thermal 
average of an infinitely small variation of the observable X. In this way, we can convert 
the Langevin equation for dynamical variables to the Fokker-Planck equation for the 
distribution functions. Here we show an explicit derivation of the Fokker-Planck equation 
for the simplest Langevin equation (2.1) as a demonstration. Eq.(2.1) can be considered 
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as a set of first order differential equations for two variables x and v — x. Then the 
Kramers-Moyal coefficients up to the second moments are given by 

C^{x)^v 

Ci{v) = V — 

m m ox 

C2{x) = 

C2{v) = ^m^ — J^ dtij^ dt2{v{ti)v{t2))\x(t)=x 

hm fi- r'^'dh^ + OiAt)] 
\t^o \At Jt w? J 



2-fT 



(2.7) 



Higher order coefficients vanish in the At — > limit. Now we get the Fokker-Planck 
equation corresponding to the Langevin equation (2.1); 



+ ( ) ■ (2.8) 



r / 7 I dV 

dtP{x, V, t\xQ, vq, 0) = dx {-vP) + dy [ —v-\ — 

[ \m m ox J ' \m^ 

This Fokker-Planck equation has a simple solution 

P«t oc e"^(^™'''+^^'^^). (2.9) 

Note that both of -vd^P + ^^d^P and [^vP + ^9„P] cancel for P'\ It is the well- 
known Maxwell- Boltzmann distribution for a system in an equilibrium with temperature 
T, and satisfies the stationarity condition dfP^^ — 0. The solution satisfies the equihbrium 
condition, stronger than the stationarity condition. 

Here we have used the words "stationary" and "equilibrium" in the following sense. 
Stationary distributions are solutions to the Fokker-Planck equation satisfying dtP = 0. 
Equilibrium distributions are also stationary but satisfy a stronger condition which is 
called the detailed balance condition. The most direct definition of the detailed balance 
condition is given in the language of the Master equation. Due to the definition of sta- 
tionarity, P'^ satisfies / dX' [w{X' X)P'\X') - w{X X')P'*(X)] = for arbitrary 
X. On the other hand, the detailed balance condition is defined as 

yx, X', w{x' x)p'\x') - w{x x')p'\x) = 0. (2.10) 

To satisfy this condition, the system must have the microscopic time reversal symmetry 
and can not have a specific arrow of time. In other words, there is no entropy produc- 
tion. In a stationary but non-equilibrium configuration, there is a fiow of current in a 
configuration space {x,v). 
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The solution of the Fokker-Planck equation can be represented in a path integral form 



as 



P{x,t\xo,0) = 



x{t)=x 



Vxexp 



x{0)=xo 



(2.11) 



Its derivation is explained in the appendix A. The "Lagrangian" L — -^{mx + ^x + ^y 



is called the Onsager-Machlup function [14]. A variation of the Onsager-Machlup function 
gives the most probable path in the stochastic processes. Apparently, since we have L >0, 
the paths satisfying L — are most favored if exist. 

The Onsager-Machlup function can be divided into two parts, 



imx+^)^ + ^x^ 



(2.12) 



which preserves time reversal symmetry, and a violating term. 



(2.13) 



The latter plays an important role to prove the fluctuation theorem in the next section. 



3 Non-equilibrium Identities 

The stochastic equations such as the Langevin or the Fokker-Planck equations describe 
how a system is dynamically relaxed to a stationary or an equilibrium state. Furthermore 
we can calculate transition amplitudes of a system to one state to another. By using the 
method reviewed in the previous section, we can calculate a ratio of an entropy decreasing 
probability to an entropy increasing probability. Since the latter probabilities have always 
much bigger values, the entropy is always increasing after wc take a stochastic average. 

In this section we review a derivation of the fluctuation theorem and the Jarzynski 
equality from the stochastic equations. 



3.1 The Fluctuation Theorem 

The fluctuation theorem was first discovered in a numerical simulation [9] and gives the 
ratio of probabilities of an entropy increasing process to that of a decreasing one. The proof 
of the fluctuation theorem is given for various systems including classical Hamiltonian 
dynamics [15], stochastic Langevin dynamics [16] and quantum mechanical evolutions 
[17, 18]. The Jarzynski equality [10] is a relation between non-equilibrium work and 
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equilibrium free energy difference, and both of them are remarkable discoveries in the 
recent developments of non-equilibrium statistical physics. In this paper, we concentrate 
on a system that the evolution is described by a Fokker-Planck equation such as eq.(2.8). 
The fluctuation theorems can be simply derived and the meaning of entropy production 
(or a violation of time-reversal symmetry) is clear. 

We consider a stochastic system described by the Langevin equation (2.1) or the 
Fokker-Planck equation (2.8). In order to study a dynamical evolution, we introduce an 
externally controlled parameter Af in the potential V{x;X[). By changing the external 
parameter Af as a function of t, the corresponding stable state changes accordingly with 
time. For later convenience, we call the process of changing the external parameter with 
Af as the "forward protocol". For example, we may set the minimum position of a 
harmonic potential as the externally controlled parameter; 



if the position moves linearly in time t, the parameter is given by Af = Vot. We can also 
take different protocols e.g. oscillatory or pulse-hke etc. 

From the path integral representation of the transition rate (2.11), a probability that a 
sequence of conflgurations = {x{t),t e [0, t]|x(0) = Xini, x{t) — xgn} is realized during 
the time interval t & [0, r] is given by 



The trajectory F,- represents a sequence of configurations in the forward protocol Af with 
the initial configuration x{0) = Xini. 

We now define a time reversal of the forward protocol Af , and call it a "reversed 
protocol" Af = Af_j. We consider a probability P^[r* |xfin] that the system experiences a 
reversed trajectory F* = {x*{t) = x{T — t),t G [0, r] |a;*(0) = Xfin, x*{t) = x^ai} in the time- 
reversed protocol Af . The reversed trajectory has the initial value x*(0) = xgn = x{t), 
x*{Qi) — —x{t). If the system has time- reversal symmetry, the probability should be 
the same as the probability P^[FT-|xini]. But since the stochastic equation violates the 
symmetry, they will be different. The reversed propability P-'^[F*|2;fin] is similarly given 



V{x-\!) = -Kx-\!)\ 



(3.1) 




(3.2) 



by 




(3.3) 
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In the last line, we change a variable from t to t' — r — t. This change causes a flip of the 
sign of X. The ratio of and now becomes 

P [rrl'^ini 



exp 



dtx (mx + ^^-^ 



dx ) 



(3.4) 



This gives a key property to prove the fluctuation theorem. Time-reversal symmetric 
terms are canceled between P^ and P^, and the ratio is given by the entropy production 
S of the stochastic process. 

We further need to sum over the initial configurations, Xini and a:fin respectively for 
the forward and the reversed protocols, with appropriate statistical weights. Here we 
assume that the external parameter is kept fixed at the initial value of each protocol 
before t = 0. Hence the system is in the equilibrium. We therefore multiply P^ or P^ by 
the Boltzmann weight P'^'^(a;ini) or ^^'^(xfin). The ratio of the Boltzmann weights for the 
initial configurations is given by 



P^q(g:ini) _ Z(Af ) 

P-l(Xfin) Z{X^) 



exp 



exp 



T V 2 



dt mxx + X 



:M^mi - %n) + y{Xim; Aq ) - V {x&n] 



^dV{x;X(}^^^dVix;\f 



dx 



dX[ 



AF' 



(3.5) 



where AF is a difference of the free energies P(A) = — T'logZ(A) of equilibrium states at 
A = A^ and A = Af , 



AP = P(Af)-P(Ao^). 
Combining the two ratios eq.(3.4) and eq.(3.5), we get the following relation, 

P [r^lXjjiijP ^(Xjni 



p^[r*\x^^]p-^{x^^) 

Here we have defined P[r,-] and VFfr,-] as 



p[r. 



1 

T 



dtx 



pdV{x; X 



dXf 



exp(P[r,]). 



(3.6) 



(3.7) 



T 



T 



(3.8) 



which measures the entropy production in the trajectory and the work exerted on the 
system. 

As a simple example, for the potential V{x; Af ) — k{x — vot)'^/2, we have 



^[rr] = -f [ dtVQk{x{t) - Vot). 



(3.9) 
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X 



Figure 1: (a) A schematic illustration of motion of a particle in a potential V{x; Af ) = 
^k{x — vqIY . This picture shows a natural configuration with {x{t) — v^t) < 0. It gives a 
positive value of -R[rT-]. (b) A noise ^ rarely pushes a particle to the opposite side beyond 
the minimum point x{t) = v^t. Since {x{t) — v^t) > 0, it gives a negative value of i?[r^] 



The term, velocity times force, gives a work exerted on the system. If we neglected the 
fluctuation of the particle, x{t) — v^t would always have a negative sign, and i^fr,-] would 
always increase. It is consistent with a naive picture. However in a mesoscopic system, 
fluctuations can grow larger and x{t) — vot can have a positive sign. Then the particle 
overshoots the equilibrium point dxV = to the positive side and i?[r^] becomes negative. 
Such a negative value of Rlr^] indicates that the system exerts work onto outside and it 
gives a negative entropy production. 

Prom the equation (3.7), by integrating all the paths of the configurations, we can 
derive the fiuctuation theorem in the final form as 

p^{Rr)= J Pa:P^[r,|Xini]P^'^(a:ini)5(i?r-i?[rr]) 

= e^^ j PxP«[r* |xfi„]p^^(xfi„)5(p, + R[r*]) 

= e^^p^(-P^). (3.10) 

The first fine is the definition of p^{Rr), i.e. the probability to get the entropy production 
Rr within the interval [0, r]. We use the relation (3.7) in the second fine. In the third 
equality the relation P[r*] = — P[r^] is used. Since the quantity Rr measures the entropy 
production in the interval, we see that entropy decreasing probabilities are related to 
increasing ones. They are exponentially suppressed, but exist with nonzero probabilities. 
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3.2 The Jarzynski Equality 

By integrating the fluctuation theorem over the entropy production, we can construct an 
equahty, so called the Jarzynski equality [19]. 

/oo /*oo 
dRrP^{Rr)e-^^ = / dRrP^{-Rr) 
-oo J —oo 

^ (e-«-) = 1. (3.11) 

We have defined the average as 

/oo r 
dRrP^{Rr)F{Rr) = / VxP"" [T r\x,^^P'%X,^0 F {R[r r]) ■ (3.12) 
-oo J 

The Jarzynski equality (3.11) states that the weighted sum of over all possible non- 
equilibrium processes with an externally controlled potential gives an unity. In terms of 
the work exerted on the system VF[r^] and the free energy difference, we can relate an 
average work done in non-equilibrium processes to the equilibrium free energy difference 
[10] as 

W AF 

(e-^) = e-^. (3.13) 

From this, by using the Jensen inequality (e^) > e^^\ we get an inequality; 

(W) - AF > 0. (3.14) 

This indicates the second law of thermodynamics. The Jarzynski equality simply states 
that there must exist microscopic processes with a large negative entropy production to 
satisfy the equality, and the probability is characterized by the equilibrium quantity of 
the free energy difference. 

Some comments are in order. First the notion of entropy is usually dcflncd for a ther- 
mal system after taking an average. So it may be appropriate to use a word, an entropy 
function, instead of the entropy for each microscopic conflguration. The second comment 
is that in the above derivation of the second law we have implicitly in mind that the above 
free energy difference is the difference between the initial and the final free energies. It is 
justified if the system is relaxed to an equilibrium state with the external parameter at 
t — T after a long time interval. Since the system is in contact with a large heat bath with 
temperature T, the relaxed state coincides with the equilibrium state at the temperature. 
If this is the case, the second law of thermodynamics is derived from the Jarzynski equal- 
ity. In the present proof of the fiuctuation theorem, we have used the stochastic approach 
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and the system explicitly violates the time-reversal symmetry. Then such a relaxation 
can occur. But if we starts from the original unitary quantum mechanical evolution, the 
system cannot be thermalized in an exact sense. In applying the fluctuation theorem to 
the information paradox of a black hole, such considerations are inevitable. 

An alternative expression of the fluctuation theorem is obtained by using a generating 
function. We define the generating function for R^- as 

Z'^iar) = In dRre"'^'^-p^{Rr)^ . (3.15) 

Derivatives of Z^{ar) give connected correlators of the entropy production i?,- in a situ- 
ation of the forwardly varying parameter. One easily gets the following relation between 
Z^ia^) and Z^ia-r) from the fluctuation theorem as 



Z' {ar) = In (ii?^e'"-«-e^-p«(-i?, 

(POD N 
J (ia;e'"('-"-V^(a^) 



^Z'^ii-ar). (3.16) 

We have used the equation (3.10) in the first line. In the second line, we changed a variable 



Rr to X = —Rr- If the forward and the reversed protocols are identical i.e. Af = A,^ 
we get a simpler relation Z{ar) = Z{i — ar). 

Finally, we give a comment on our assumption of the initial distribution. We have 
assumed that the initial distribution is an equilibrium one. This condition can be easily 
relaxed to a steady state. More generally, if the initial distributions for Xjni and xgn arc 
-P^n^ini) ^-iid P^^{xfin) respectively, we can define an entropy production as 

R[^r] ^ In ( C^NS^l • (3-17) 



T-tl 



VP«[r;|Xfin]P«*(^fin) 



Then we get the fluctuation theorem in the form; {R^) / p^{—Rt) — e^'^ . The choice of 
initial distributions is arbitrary, but the problem is that we usually do not know an explicit 
form of the distribution function of a steady state P^*. The steady state fluctuation 
theorems are reviewed in the appendix C. By using it, we can derive the Green-Kubo 
formula and its non-linear generalizations. 
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4 Langevin equation in a Black Hole Background 



In this section we derive a stochastic equation for a scalar field in the black hole back- 
ground. We take /i — )■ hmit with the Hawking temperature Hk /2t^ fixed. Since the 
energy is absorbed into the black hole, a dissipation term is induced at the horizon. The 
classical equation is furthermore modified by the quantum effect, i.e. the Hawking radi- 
ation from the black hole. Because of these effects, the equation of motion in the black 
hole background is described by a stochastic Langevin equation with a quantum noise 
and a classical dissipation terms. We first review the basics of black holes and the Hawk- 
ing radiation, and then derive the Langevin equation of a scalar field in the black hole 
background. 



4.1 Space-time Structure of Black Holes 

First we summarize some basic facts of the space-time structure of black holes. (For a 
review, see for example [20].) Here, we consider a spherically symmetric neutral black 
hole, the Schwarzschild black hole. It is a solution to the Einstein equation in vacuum 
with a zero cosmological constant and the metric is given by 

fir) 

fir) = 1 - dQ^ = dd^ + sin^ dd(t)''. (4.1) 

M is the mass of the black hole and the only parameter of the solution. The solution is 
asymptotically flat; it approaches the flat metric at the space- like infinity r — > oo. It has 
time-translation symmetry and the associated time-like Killing vector is given by ^ — df. 
A Killing horizon is defined as a null hypersurface on which there is a null Killing vector. 
In the present case, it is given by the condition g{^,^) — —f{r) — -ir^ r — vh — 2GM. 
The surface gravity k is defined on the Killing horizon via the relation 

= (4.2) 

A direct calculation shows that k, = f'{r)/2\r=rf, = 1/4G'M for the Schwarzschild black 
hole. 

There are several different definitions of horizons. An apparent horizon is a more 
general concept and defined locally as the most outer trapped null surface. It does not 
need a time-like Killing vector as the Killing horizon, but it is defined in an observer- 
dependent way. An event horizon is defined in a global way as a boundary of the past 
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light cone of the future infinity. MathematicaUy a black hole is defined as a set that is not 
contained in the past light cone of the future infinity. For the Schwarzschild black hole, 
all the definitions of the horizon coincide, but they are different for dynamical black holes. 
In applying non-equilibrium statistical physics to the dynamics of black holes, we need to 
pay special attentions to the differences. In the present paper, however, since we consider 
an eternal black hole as a background space-time, their differences are not essential. 

The coordinates used in eq.(4.1) is called the Schwarzschild coordinates. The singu- 
larity of the metric at the horizon r — rn is not physical, and can be removed by using 
other coordinates, such as the Kruskal (-Szekeres) coordinates (C/, V) 

U = --e-'^it-r,) ^ V = ie''(*+"*) (4.3) 

K K 

/dv T 
-— =r + rj/log( 1). (4.4) 
f{r) th 

Here r* is the tortoise coordinate and takes — oo < r* < oo between the horizon and the 
spacial infinity. In terms of the Kruskal coordinates, the metric of the Schwarzschild black 
hole becomes regular at the horizon; 

ds'^ = -—e'^dUdV + r^dQ^ (4.5) 
r 

At the price of removing the coordinate singularity, the asymptotically fiatness is unclear 
in these coordinates. We will impose regularity conditions on physical quantities at the 
horizon in the Kruskal coordinates. 

Figure 2 is the Penrose diagram of the Schwarzschild black hole, which captures the 
causal structure of the space-time. The vertical and horizontal axises correspond to the 
Kruskal time T ^ {V + U)/2, and the Kurskal radius R ^ {V - U)/2. In contrast to 
the Schwarzschild coordinates, the Kruskal coordinates are regular beyond the horizon 
(r = rn), and can be extended to the maximally extended Schwarzschild space-time 
(— oo < U,V < oo). The original Schwarzschild coordinates (— oo < t < oo^rn < r < 
oo), on the contrary, can cover only the region I in fig.2. We define (i, r*) coordinates 
in other regions. For example, in the region n, we can define them by the relations 
U — e~'^^*~^*y K,, V — —e'^^^~^'^*y K,. In the Kruskal coordinates, the space-time is separated 
by the future and past event horizons {U — and V — respectively) into four regions. 
There are four possible combinations of signature of U and V as shown in the table 4.1. 

Finally we note that the time-like Killing vector ^ = dt is written as ^ = K{Vdv — 
U du) — nRdr in the Kruskal coordinates and, therefore, the directions of time are opposite 
in the region I and II. We have drawn the directions of ^ in fig.2. 
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U r = V 




Figure 2: A point in the diagram represents a two dimensional sphere with radius r 
at time t. r-constant and i-constant surfaces are depicted. Arrows on the r-constant 
surfaces indicate the flow of the time-hke KiUing vector. They have opposite directions in 
the region I and II. The singularity at r = is drawn by zigzag lines in the diagram. Event 
horizons are located aX, r — th and separate the space-time into four distinct regions, i"*", 
i~ and Iq are the future, past and spatial infinities. 



I 


[/ < o,y > 


r>rH 


n 


[/ > o,y < 


r>rH 


in 


[/ > o,y > 


r <rH 


IV 


[/ < o,y < 


r <rH 



Table 1: Four regions of maximally extended Schwarzschild space-time 



4.2 Field Theory in the Black Hole Background and the Hawk- 
ing Radiation 

We briefly summarize the quantum field theories in the black hole background. For a 
comprehensive review, see e.g. [21]. The action of a massive scalar field in the maximally 
extended Schwarzschild space-time is given by a sum of the fields in the right wedge 
(region I) and in the left wedge (region IE). In each region, the action is given by 

(4.6) 
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where we have decomposed the field into partial waves 

0(t,r,Q) = 5]^^^^^n„(Q), (4.7) 

l,m 

and defined the effective potential for each partial wave with an angular momentum Z, 

V;W = /w(^ + ^^ + m^). (4.8) 
The equation of motion of the scalar field is given by 

[dl-dl+Vi{T)\4>RMl,m)^^. (4.9) 

Both in the asymptotically flat region (r — )■ oo) and in the near horizon region (r — )■ r^^), 
the potential Vi vanishes and the equation of motion is reduced to the free field equation. 
Thus, in the near horizon region, the classical solutions are approximately given by 

{ 1 ^-iixiut+ikr* T>\ 
^ ,1 (4.10) 
(m L) 

T, . \ (inR) , ^ 

' ^ I _i=e*^fc*+*^''* (InL). ^ ^ 

and their complex conjugates. Here Uk = +\k\ > 0. The sign difference in front of iojkt 
in R,L follows the convention of [21]. With this convention, these fields are positive 
frequency modes with respect to the time-like Killing vector, dt in R and —dt in L, 
satisfying C±d^Uk = —icOkUk- The complex conjugates {u^'^)* are the negative frequency 
modes (in the above sense) satisfying C±q^uI. = +iUkU^- They are orthonormal with 
respect to the following Klein-Gordon inner product. 



(/, g)^t [ d?x^, {fdtg - dj*g) 

j dr^ {f*l,rn)9t9{l,m) - dtf*i,m)9{i,m)) ■ (4.12) 

l,m 

The integration is performed on a constant time slice E^, but it can be generalized to any 
space-like surface E and the choice of the integration surface does not change the value 
of the inner product. The field (f)(i^rn) can be expanded in terms of the classical solutions 
in the Schwarzschild coordinates in the near horizon region as follows; 

/dk 
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We will omit the suffixes (/, m) of creation and annihilation operators for simplicity in the 
foUowings. 

In the Kruskal coordinates near the horizon, the equation of motion becomes dudv4>(i,m) 
{d^ — 9^)0(/,m) = 0. So we may define another basis of functions 

<(r, R) = (4.14) 

where Ep — +\p\ > 0. They are positive frequency modes with respect to the Kruskal 
time. In terms of them, the field can be expanded as 

h,m) = J M + Mi<r] ■ (4-15) 

In contrast to the wave functions (4.11), they are defined globally in the whole space-time. 

In order to relate two different definitions of the creation and annihilation operators in 
the Kruskal and Schwarzschild coordinates and to express the Kruskal vacuum bk\0)K = 
as a Fock state constructed on the Schwarzschild vacuum a^'^\0)R,L = 0, we look at 
the analyticity properties of the functions [22]. The positive frequency wave function 
Up in the Kruskal coordinates with p > (or p < 0) is an analytic function in the 
lower half U (or V) plane since ~ ^-iEpU ^ ^-lE^vy -^.j^g other hand, 

since e*'^^* = {r/rn — l)*^e*^'", there is a phase jump when it crosses the horizon. So we 
need to combine the positive and negative frequency wave functions in the Schwarzschild 
coordinates to construct a wave function with the same analyticity property as . They 
were obtained by Unruh [22] as 



(1) 1 f R , _:i!fk 



h 

(2) 



•^/2sinh 



These combinations are analytic in the lower half plane of C/ or V. In the following we 
set h — 1 for notational simplicity. Such analyticity property can be easily checked. For 
example, u^^'' with a positive k can be rewritten as an analytic function of U 

oc (-kC/)^, (4.17) 

if it is analytically continued from the region I of the right wedge (t/ < 0) to the region II 
of the left wedge (?7 > 0) through the lower half of the U plane by the transformation U — > 
C/e"^. Hence the combination is analytic in the lower half plane of U. For A; < 0, u'^^ oc 
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{kV)~ and it is analytic in the lower half plane of V as e~'^*^. In the classical hmit 
where /i — )■ 0, u~^^ becomes a positive frequency mode in the Schwarzschild coordinates 
g-«wfc(t=Fr,) locahzed in the region I. Similarly, u~^^ with a positive momentum > is 
written as an analytic function of the lower half plane of V, {nVy^'^f'^ while, for a negative 
k, it is analytic in the lower half plane of U and written as (— acC/)~*'^'=/'*. It behaves as 
a negative frequency mode in the Schwarzshchild coordinates but localized mostly in the 
left wedge in the classical limit. They penetrate into the right wedge by quantum effects. 
In this sense, u^^^ is classical while li^f ^ is quantum in the right wedge. 
The scalar field can be expanded in terms of these modes as 

ni,m) - J ^y^^ [^k % + (Cfe j + Cjfc % + [Ck J'K J J • (4.18) 

The Kruskal vacuum (6p|0)x = 0) is equivalcntly given by the conditions, c^j}^\0)k = 
c^^^\0)k = 0. The annihilation operators in the Schwarzschild coordinates and can 
be expressed as a linear combination of c[,^'' and c^^^ as 



y/2 sinh - 



1 f (2) 

' e 2. c) 



(4.19) 

where 71(0;^) = i/(e^'^'*'*^/'^ — 1). Hence the Kruskal and the Schwarzschild operators are 
related by the Bogoliubov transformation. 




(1) 



+ ^(^fc) V^(^fc) \ [ ('k' \ _ rr [ 



The transformation can also be represented as 

sinh^^fe = n(6^fe). (4.21) 

From this transformation law, we can read off the relation between Kruskal vacuum and 
Schwarzschild vacuum as 

|0)k = e-^^|0)fl|0)^ (4.22) 
^ 00 

A; n=0 
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Note that the Fock space in the left wedge \n^) is constructed on a Minkowski vacuum 
with the backward time direction (— t). 

The expectation value of the Schwarzschild number operators (af )^af in the Kruskal 
vacuum \0)k is given by 

x(0|(af)taf|0)x = ^— ^e-^ k{0\^1&\0)k 

K 

^n{cuk). (4.24) 



e " 



This is the thermal distribution of the Hawking radiation [2], and characterized by the 
temperature = nh/2Ti. Note that the thermal spectrum in the right wedge is created 
by the effect of the field u\. , which is classically localized in the left wedge but penetrates 
into the right quantum mechanically. 

For a generic operator 0[i = O^^a^, ('^^)^) which is made of only and {a^)\ its 
expectation value k{^\Or\Q)k can be interpreted as a thermal average. Such a thermal 
behavior can be generalized to products of operators, such as k lO k , made of 
both the right and left creation (annihilation) operators. Its expectation value can be 
interpreted as a Schwinger-Keldysh correlator. 

First let us consider k{^\Or\Q)k- Since the Kruskal vacuum is represented as (4.21), 
one has 



k ^ n=0 



= Tr 



R 



o 



R 



(4.25) 



Zr 

Here, the Hamiltonian and the partition function are defined by 
Hr^ j ^a;,(a^)taf , 

oo 

Zr = Tr [e-'^^^j = JJ e"^^^"^ = cosh' dk- (4.26) 

k n=0 k 

Hence k{^\Or\Q) k can be interpreted as a thermal average of the operator Or at the 
Hawking temperature Th- 

For a product of left and right operators, the expectation value in the Kruskal vacuum 
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is given by 

^ oo 
fe,fe' ^' m,n=0 



(4.27) 



cosh cosh Ok 

k,k' " m,n=Q 



In order to express it as an expectation value in the right wedge Fock space, we first 
rewrite the expectation value {n^y.\0 L\m^^,) in terms of the operator in the right wedge 
as follows, 

(m^,|(9^|nf:,,) = (nf,|6^|mf), (4.28) 
Here we have defined the operator O^aR, a^) by the following substitution, 

d,{a\ (a^)t) = J2'^n.M)"\4^r ^ Olia^a^^) ^ c^,n(a^,)"(«-l)"'- (4-29) 

n n 

Note that the coefficients c„j „ are not converted to its complex conjugate. In particular, 
the field itself ^^(t) is converted as 

/flk 
[a^e''"'=*+''="* + {a^)^e-"^>''-"''-*] (4.30) 

^ (PKt) - j JL[(a;Rjteia;,t+ifer, ^ ^R^^-u.,t-ikr,^_ (4 3^) 

This has the same functional form as (jy^it). We later interpret this field 0^(t) as the 
(lower) Schwinger-Keldysh field and write as 0r to distinguish the original right-wedge 
field (f)R. By this substitution, the expectation value can be written as 

K{^\dRdL\^)K = ^Tr (e-t^« Or e-i^« ofj (4.32) 
= {ORdl),, 

2 ' 2 

In order to distinguish it from the ordinary finite temperature Green function, we have 
introduced the notation {■■■)§_ n as above. 

2 ' 2 

If Ol is made of a product of operators Ol = AlBl, it is converted as 

Ol ^AlBl^OI^ Bl Al. (4.33) 
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Especially a care should be taken for the time evolution operator UL{t,to) = Texp 



i J^^ dtHL{t) ■ Following the above substitution rule, it is converted to 



UL{t,to) = Texp 



-I I dtHLit) 

'to 



Ul{t,to)^fexp 



t 

-I I dtHl{t) 

'to 

V 



(4.34) 



where T is an anti-time ordering. For a hermitian Hamiltonian, = Hr is satisfied and 



t 

-I I dtHnit) 

'to 



(4.35) 



Hence a Heisenberg operator Al{x) is mapped to 



M^) = Ul{t^,to)AL{to)UL{t^,to) ^ Al{x) = U^{t^,to)Al{to)Ul\t^,t 



(4.36) 



The converted Heisenberg operator is evolved backward in time with the Hamiltonian 
(-Hr). 

From these considerations, an expectation value of a general operator including both 
of left and right operators can be represented as a path integral form of the right-handed 
fields; 



/3 /9 



/ 



V(f>RV^R OR[(f>R] OH^r] eXp iS[<Pr] - iS[^R^ 



(4.37) 



Here (pR represents the original right-wedge field while a new field (f>R{t) is introduced to 
represent the transformed operator (9^. The minus sign in front of the action S[(j)R] comes 
from the backward time-evolution of O^. If we combine 0i?(i) and 0i?(t) together as a 
single (f)R{t) field along a doubled path depicted below, this expression is equivalent to the 
closed time path formalism of the real-time finite temperature field theory. The insertions 
of exp(— /3i?R/2) can be represented as an evolution of time into the imaginary direction 
with —(3/2 at both ends. Hence the path is given on the complex time plane as Fig. 3. The 
field on the lower line corresponds to the field in the left-wedge as (/)R{t — i/3/2) = 0l(^)- 

An alternative interpretation is an analogy with the thermo field dynamics [23], an- 
other method to deal with the real-time finite temperature field theory. In this analogy, 
the operators in the left wedge can be regarded as the "tilde-fields" of thermo field dy- 
namics [24]. 
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— oo 



Or 



OO 



i 



— oo — i/3 



oi 



oo-f 



Figure 3: lives on the upper line while (pR lives on the lower line. The time evolution 
is backward on the lower line. 

4.3 Effective Equation of Motion in the Vicinity of the Horizon 

Now we derive an effective equation of motion for the scalar field in the black hole back- 
ground. Classically a dissipation term is induced since the energy is absorbed into the 
black hole horizon. In quantizing the system, a noise term will also be induced because 
of the Hawking radiation, and the system is effectively described by a Langevin equation. 

The effect of the absorption can be described by imposing the ingoing boundary condi- 
tion at the horizon r = Th. Since, in the near horizon region, the system can be described 
by a set of 2-dim free fields satisfying (c?| — d^J(j)(^i^rn) = 0, the ingoing boundary condition 
can be represented as 



The condition implies that there are no outgoing modes at the horizon, and violates the 
time reversal symmetry. 

Since the scalar field is coupled to the gravitational field, if it is quantized, the chiral 
condition at the horizon seems to violate the general covariance by the quantum grav- 
itational anomaly. The violation is compensated by the fiux of the Hawking radiation 
[25, 26, 27]. In the following we will see that the quantization of the scalar field near 
the horizon naturally leads to the chiral condition (absorption) with the flux of Hawking 
radiation (noise term) at the horizon. 

The method we will use is similar to the retarded-advanced (or Schwinger-Keldysh) 
formalism. The derivation of a Langevin equation is given by integrating fluctuating 
fields. (For a review, see, e.g. [28].) 

4.3.1 Integrating Out the Environments 

In obtaining the Langevin equation at the horizon, we need to integrate out certain kinds 
of environmental variables interacting with the system variable at the horizon. In order to 
do this, we first consider a stretched horizon at r = r^/ + e and treat the variables between 
the horizon (r = rn) and the stretched horizon (r = ri/+e) as the environmental variables. 
Because of the quantum mixing of the wave functions in the left and right wedges (4.16), 



{dt - dr,)(t)[i,m){t, r^rn)^ 0. 



(4.38) 
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the integration of these variables corresponds to an integration of the fields in the left 
wedge, which are classically hidden. In this way, we derive a Langevin equation at the 
stretched horizon. This equation is shown to be independent of the small parameter e 
characterizing the position of the stretched horizon and we can take e — > limit at last. 

Since the Langevin equation we are going to derive is the equation of motion at the 
boundary of a rigion [rn, th + g], it is convenient to discretize the equation of motion near 
the horizon. In the tortoise coordinate r* in which the equation of motion becomes free, 
the region is mapped to [—00, r*{rH + e)]- We divide the region into infinite segments as 
('"*)n = r*{rH + e) +nd (for n — 0, —1, —2, • • • — 00) and set oscillators x„ on these lattice 
points. Here d is a lattice spacing in the tortoise coordinate. Discretized equations of 
motion for the scalar field are given by 

Xq = -k{xo - xi) + k{x_i - Xq) - V;((r*)o)a;o 

; (4.39) 

X—n — k{X—n X—n—l) + k{X—n-\-\ 2^— n) ^((^*)— n)^— n- 

The continuum limit is given by taking d — >■ limit with kd? — 1 and (t>{i,m){t-i ('"*)«) — 

Xn{t)/Vd. 

k k k k k 

^^^^ 



X_2 X_i Xq 



Xi 



Figure 4: Discretized model of a scalar field in the near horizon region. The variable Xq 
represents a variable at the stretched horizon. 

Introducing the forward and the backward differentials in the tortoise coordinate r*, 
A+x„ ^ ^!^±1^ , A-Xr, ^ ''--^--\ (4.40) 

and using the relation, 

-{Xn - Xn+l) + {Xn-i -Xn)^d (A+ - A") x„ = d^A+A'x^, (4.41) 
we can write the discretized equation (4.39) for n < as 



(^i,m{{n)n) = kd^A+A (l)i,m{{r*)n) " VJ((r*)„)0i,^((r*)„) 



> 0;,m(r*) = d;ji,m{n) - Vi{n)(j)i,rn{n). (4.42) 
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Vi{{r^)n) stands for the gravitational potential at (r*)„. It is proportional to /(r) = 
1 — 2M/r and vanishes in the near horizon region. Hence we neglect the potential term 
later in this section. 

The normalization of the fields 0(i,m)((^*)n) = Xn/Vd is determined from the action. 
With this normalization, the discretized action becomes the continuum one; 



1 " 

L n=— oo 



- / dt 



dr^ 



(0a,m)(r.))'-(5,.0(,,^)(r,))'-V;(n)(0a,„z)(n))' • (4.43) 



Here the discretized potential U is defined by 



(4.44) 



Note that, we define dJl ^ ^ r*(xH+e) ^ ^ ^ ^ 

' Z— /n=— oo J— oo 

The full action is a sum of the fields in the left and the right wedges. As we saw 
in the previous section, the path integral containing both the left and right fields can 
be rewritten by a path integral of a right field on a closed time path. In the previous 
section, we have written the filed in the lower line by xr. In the following, we use a unified 
notation and write xr by x^ and xr by x^. 

Previously we considered a path from t — — oo to oo. It can be generalized to a path 
up to a finite time with fixed boundary conditions xl{t) — Xf^^ {1 — 1,2); 



'[4n,^]= / 



Y\. ^'^^^'^^e''^t^-^'"''^o'^il~*'^'^-JV"'^o.^il (4.45) 



-oo 



xt 



't - 



if} 



-oo — IL 



Figure 5: The values of the fields at the right ends of the paths, a;^(t) / = 1, 2, are fixed 
in the path integral. 

By integrating the environmental variables (fields between rn and th + e), we have 

" P2;^P2;2e^^[^o'^5I-'^[^o'^i]+*^^^[^o'^ol. (4.46) 
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The definition of the infiuence functional Sjp is schematically written by 



e 



where Sint[x^i,XQ] = k/2 f dt{xQ — a;_i)^. 

Since the system variables Xqjx"^ are coupled linearly with the environment variables, 
the infiuence functional Sjf[xI,xI] has a Gaussian form 

Sif[xI,xI] = dtdt'xl,{t)Fij{t,t')xi{t'). (4.48) 

The Kernel function in the Schwinger-Keldysh formalism Fjj{t,t') can be obtained by 
taking derivatives of the infiuence functional as 



\xl=0 



1 P 

Ft lit f') — - - piSiFlxi,xl 

' > ~ i 8xi{t')8xi{t) 

^ / {TA^xUt)A^^Ut'))ii -(A+xi,(t)A+a;^(t'))|, 

^ [ -{A+xl,{t)A+xl,{t')), , {fA+x\{t)A+xl,{t')), , ^ - ^ 

\ 2 ' 2 2 ' 2 

The expectation means an integration over the environmental variables x^_^,--- 
(/ = 1,2). T stands for the time ordering, and T is the anti-time ordering. As we saw 
in the previous subsection (4.37), these propagators are equal to the Schwinger-Keldysh 
(SK) ones with the path drawn in Fig.3. In the continuum limit, the discrete Green 
functions d x Fij{t,t') become the continuum counterpart 

-^r.Cri I 21 (t r -t' r') (t r 't' r') / ^'='' -'"+' ^ 

and the influence functional is given by 

SiF[<i>\rH + e), 0'(r^f + e)] = ^ J dtdt'c^l^^it, rn + e)Fli^,mUV , rH + e). 

(4.51) 

Strictly speaking, the expectation in the Green functions should be evaluated at (r*)_i, 
but in the continuum limit it coincides with the position at the stretched horizon at th+g- 



\r=r'=rH+e 
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4.3.2 Vacuum Condition 

From the previous discussions, we already knew that the Green functions in the Kruskal 
vacuum become identical with the Schwinger-Keldysh Green functions along the contour 
in Fig.3. We repeat the discussion for the case of the two point functions explicitly in the 
following. In order to calculate the influence functional, we need to specify the vacuum 
condition for the environmental variables, i.e. the Kruskal vacuum condition so that the 
physical quantities is regular in the Kruskal coordinates. We expand the scalar field by 
u^^\u'^^ and its complex conjugates 

Sn .(t r)- f \c^^^ M^'^ + fc^'^ )Uu^^^Y + c^^^ M^'^ + fc^'^ )Uu^^^Y 

(P{l,m){l^,r*) — J ^ 4^^^ [^k{l,m)'^k ^ y^k{l,m)) ) ^ ^k{l,m)'^k ^ \^k{l,m)) \^k ) 



(4.52) 



with the canonical commutation relations 



c 



(1) /Ji) 



kii,n.y i^'ii'M^f] = {'27r)2u,5u'5rr,m'Sik - k'), (4.53) 
^^m,m)^ (cg)('' -'V] = (27r)2a;,5„,5w5(fc - k'). (4.54) 

The correlators in the Kruskal vacuum become the following forms, 

FK!(i,m)ii',m')it, t') = 5ir6mm'drA'Gjf{t, r,; t\ r'X,=r', (4.55) 

where 

Here K means the expectation value in the Kruskal vacuum. As we saw, they are related to 
the Schwinger-Keldysh Green functions F^^-^ (t, t') = 6ii'6mm'9r^dr'^G''-^ {t, t') l^^^/^^^+g 
discussed in the previous section as 

^ j dtdt'(t>l^^{t,rH + e)F^%,rnUl',m')M^^^ 

= i y dtdt'cPli^^^it, th + 6)4^i),(,,,^,)(t, t')4',m')(^'' ^H + e), (4.57) 

where the Schwinger-Keldysh Green functions G^^ are given by 

" (T0i(t,r,)0i(t',r:))| | 




2 ' 2 2 ' 2 

dk 1 



47ra;fc2sinh(7ra;fe/K) \ M'^\t,t') M'^^{t,t') 



58) 
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Non-diagonal entries have an extra minus sign with respect to eq.(4.56), since the (fP field 
is defined to evolve backward in time as in Fig.3. 
Each component can be calculated as 

M^\t, t') = eit - t') (e^e--(*-*') + e-^e-(*-*')) 

+ e(t' - t) (e-^e-^'^(*-*') + e^e^'^(*-*')) , (4.59) 
M^\t,t') = eit' - t) (e^e-^-(*-*') + e-^e^-(*-*')) 

+ d{t - t') (e-^e--(*-*') + e^e-(*-*')) , 

It can be also rewritten in the following form, 

IJ(^ ^ .J J \ _ f ^^hdk ikoit-t')+ik(r,-ri) 



-k'^+ujl-ie + 27rm(a;jfc)5(-A;g + ujD -2m^/n{l + n)5{-k^ + ujI) 

-2my/n(l + n)5{-k^ + ul) z^^^J^ + 2mn(uk)5{-eo I' 

' fc ' ' 

(4.61) 

where n{uk) = l/{e^"^^ — 1), Ph = Stt/k. The 2-2 component of the Green function 
coincides with the anti-time ordered finite temperature Green function, while the 1-1 
component is the ordinary time ordered one. 

In the conventional real-time finite-temperature field theory, the contour is usually 
taken as in the figure 6. The contour corresponds to considering an ordinary finite tem- 

-o o , ^ oo 

— oo — 

Figure 6: This path corresponds to the propagators which have non-diagonal entries 
(4.65). 

perature Green function with the Boltzmann factor e~^^ at the left-end; 

(O^O')^ = -^Tr {e-^''d\t) d^{t')) (4.62) 

irrespective of whether these operators live on the upper or lower lines. On the other hand, 
the Green function (4.37) we are considering corresponds to taking a different contour as 
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drawn in Fig. 3. The fields on the lower line in these two contours are related by the 
similarity transformation 



<l>L^{t, x) = e-P^-/^\t, x)e^^-l^ = + i/3/2). (4.63) 

Here 0^ and 0ncw fields appearing in the formalism of Fig 3 and Fig 6 respectively. 
In the momentum representation, it is 

</'NewW = e^0^(fc), (4.64) 

In terms of the new field, the Green function can be written in the following form. 



{T<P\t)<j>\t')), -(0^(t)0Lw(t'))/3 



^^0^^ -iko(t-t')+ik(r,-r',) 

-2m sgn(A;o)(l + n{kQ))5{-kl + ul) + 2Tiin{uJk)5{-kl + ul) 

4.3.3 Langevin equation at Stretched Horizon 



(4.65) 



The effective equation of motion at the stretched horizon can be obtained by taking a 
variation of the effective action S[xl] — S[xq] + Sjf[xq,Xq]. In taking a continuum limit, 
a care should be taken since we have already integrated out the environmental field x^i, 
and only the interaction with the outer variable Xi appears in the effective action for xq. 
The equation of motion for Xq becomes 

= -k{x^Q -xi)- J dt'F^\t, t')xi {t'). (4.66) 

In the continuum limit {d — > 0) with kd? =fixed, the time derivative term drops and we 
have 

(Note that the discretized d x F^'^ becomes the continuum -^(;'^)(;/ ^/)-) The dynamics 
seems to have disappeared in the effective equation at the stretched horizon, but we will 
see that another time derivative term (which is first order) is induced from the second 
term. 
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In order to show this, following the retarded-advanced formalism discussed below, we 
recombine the Schwinger-Keldysh fields, (f)^{x) and into a cZassica/ variable (p''{x) 

and a fluctuating variable (f)"'{x). The interpretation of classical and fluctuating variables 
comes from an observation that the action — S[(p] has a dominant contribution in 
the path integral when the configuration of two fields coincide. As we saw in Fig. 3, 
the time axis of (j)^{t) differs from that of (f)'^{t) by an amount of /3/2 into the imaginary 
direction, and dominant configurations are given in terms of the redefined field (4.64) 
^New(0 — <f>^it + ^/5/2) in the following way. We define the classical and fluctuating fields 
as 



(4.68) 



Propagators are transformed in this basis as 



where we have defined 

G\t) = i(Gj,L - G],l^ + G'^\^ - Gi^Jit) = te{t){m, 0(0)]) (4.70) 

G^{t) = ^(GJ^L + G'J.L - G'^\^ - G^L) W = -iOi-mm 0(0)]) (4.71) 
G^y-{t) = -^(GJ,^ - G],l^ - Gl\^ + G'^^Jit) = \{{<p{t), 0(0)}) (4.72) 

and used the relation Gjj^^ + G^ew + ^New + ^New — 0- Because of this, the basis {(f)'^ (ff) 
are often called the retarded-advanced basis. 

In terms of the r, a-fields the influence functional can be written as 



dtdt' 



+#a,-)W^-^^i^;^::^)a',mo(^'^')0a',„.')(^')]- (4.73) 

The derivative of the retarded Green function d^^d^'fi^ satisfies 

drA',Gfi^^)(j^, ^^,){t, t')\r=ri=rH+t = -Sll'Smm'dt'6{t - t'), (4.74) 

On the other hand, the symmetric Green function can be written as 

C^ym = J JL. (n + (^e-*"*(*-*') + e+^'^^(*-*')) g^M-^-O^ (4.75) 
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and its derivative becomes 

dr^dr'^G^n ^n, ,Jt,t')\r=r'=rH+e — ^ll'^v 



^^0 ^0 ^-iko(t-t') 

47r tanh ^ 



U'fJmm' 



27r 



K 



4sinh 



2 K{t-t') 



(4.76) 



for t ^ t' . The integral is divergent t — t' . Since we are interested in the finite tem- 
perature effect, we regularize the symmetrized correlator by removing the ^-independent 
divergence (note T — k/2t^) as 



77T + 



4sinh2^ {t-t'Y 



(4.77) 



Hence the action for the stretched horizon variable, which is a sum of Sif and the inter- 
action term with the neighboring variable xi, becomes 



j dtdt'dQrl 



+ int,re,n)K{t,t')<j)%t',re,n) 



(4.78) 



Here, = th + ^ appears with rewriting ^^^^^ to 0^. By integrating the fluctuating 
variable 0"(i), (4.45) is written as 



t] 



T>(f)^ exp 



dtdt'dnr^dt - dr*)4)''{t, r„ n)K-\t, t'){dt> - dr'*)(f)''{t' , r„ Q) 



(4.79) 



It describes the effective dynamics at the stretched horizon. Note that the effective action 
contains a term which is odd under the time reversal transformation. 

Instead of integrating out the fluctuating variable, we can introduce an auxiliary fleld 
m by 



exp (^-J dtdt'ril,rn)mMr^l,rn){t^^^ 



Viexp ( Z I dtril,m){t)V2^(l,m){t) - ^ I dtdtXil,m){t)K-'{t,t%i,m){t') ) • (4.80) 



Then the probability to see (p' (t) = at the stretched horizon is written in terms of the 
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scalar fields and the auxifiary field ^ as 




The variation with respect to 0" gives the equation of motion for (ff 



(4.83) 



with the Gaussian noise correlation 



{i{l,m)it)) = , {^{l,m){t)^(l',m')it')) = h'^ram' K{t,t'). 



(4.84) 



As expected, if we take the statistical average, the outgoing modes vanish in the averaged 
sense {{dt — drS)(t>^) — 0, which means that there are only ingoing modes at the (stretched) 
horizon. The noise term can be considered as the effect of the Hawking radiation. In the 
next subsection, we compare the noise correlation obtained here with the fiux of the 
Hawking radiation. The noise correlation is not white, and the memory effect remains 
with a time scale of the Hawking temperature {t — t') 1/k — h/27rTH. If we look at 
the dynamics of a time scale larger than it, we can approximate the noise as the following 
white noise 



The above effective action is obtained previously based on the physical picture of the 
Hawking radiation [5] or a technique to reproduce the Schwinger-Keldysh formalism [7] 
in a setting of vibrating string in AdS black hole background. We have reproduced the 
same effective action by explicitly integrating the environmental variables between the 
horizon and the stretched horizon. Because of the mixing of the wave functions (4.16), 
the integration corresponds to an integration over the variables hidden in the horizon. 



5 The Fluctuation Theorem for Black Holes and Mat- 



Now we apply the fiuctuation theorem to the scalar field in the black hole background. 
Most generally, we must treat the whole system of the scalar field and the space-time 




6{t - t') = 6ii'6m,m'TH6{t " t') 



(4.85) 



ters 
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as a coupled quantum system, and backreactions to the space-time structure must be 
included. In our previous letter [29], we briefly sketched how to treat the metric degrees 
of freedom quantum mechanically in the path integral formalism and discussed the effect 
of the backreaction. In the present paper, in order to give a more systematic and complete 
investigation, we consider an easier situation, i.e. a scalar system in a fixed black hole 
background. We neglect effects of backreactions. Even if we adopt such a simphfication, 
various interesting results follow the fiuctuation theorem applied to our system, such as 
a proof of the generalized second law and a derivation of the Green-Kubo formula. 

5.1 Discretized Equations outside the Stretched Horizon 

The equation of motion of the scalar field (jf{t^r^) in the black hole background consists 
of the two coupled equations, namely, the effective equation at the stretched horizon 
r = r// + e and the bulk equation of motion outside the stretched horizon. 

We put the scalar field in a box with a radius rB(> th) and impose a boundary 
condition at the outer boundary (j){t,r = tb,^) = in this subsection. Owing to the 
boundary condition, the scalar field is shown to be thermalized. Another merit of confining 
the system in a box is to stabilize the total system (even we take the backreaction into 
account [30] if the size of the box is not so large.) It thus justifies to choose an equilibrium 
distribution as an initial distribution for the matter field as we will do in the following. 
In a later section, we choose a different boundary condition to realize a steady state with 
a constant energy ffux. 

In order to apply the fluctuation theorem reviewed in section 3, we need to construct a 
Fokker-Planck equation which is local in time. In doing so, it is necessary to approximate 
the noise correlation (4.84) by the white noise (4.85). This approximation is valid for 
a longer time scale than H/^t^Th. Though the validity is hmited, we consider such an 
approximation in the present paper. The memory effect of the colored noise will be 
discussed later. 

In the white noise approximation, the discretized equations are given by 

' 7oXo = -K^o - x^) - V2^o , iUtMt')) = ^oTnSit - t') 
xi = -k{xi - X2) + k{xo - xi) - Vi{ri)xi 

< : (5.1) 

xn = -k{xN - xn+i) + k{xN-i - xn) - Vi{rN)xN 
^ xn+1 = 0. 

The flrst fine is the stochastic equation for the fleld at the inner boundary (stretched 
horizon) r — = rn+e with a noise term ^. Note that the time derivative term originates 
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Figure 7: A schematic illustration of eq.(5.1). Each variables Xi are bound by spring with 
constant k. Only xq reserves a friction 70 and a noise ^o- The effects of Vi{ri) are not 
described on picture, it affects each variables as harmonic potential. 



in the dissipation term induced by the interaction with the environmental variables. The 
last one is the boundary condition at the outer boundary r — tb- The middle ones are 
the bulk equations of motion, and the field 0i,m('"*) between the stretched horizon and 
the outer boundary is discretized into N lattice points in the tortoise coordinate r*. The 
continuum limit can be taken as before by taking d — > (or N ^ 00) with the following 
conditions and replacements; 

{N + l)d^r^{rB) -r^{r^), kd^ = 1 , 70^= 1, 0(i,^)((r*)i) =Xi/y/d, {(;,^) = ^{0, 

(5.2) 

where d is the lattice spacing in the tortoise coordinate r*. The continuum equations in 
the bulk can be recovered as before and becomes 

^(/,m)(^> r*) = dl,(l>\i,m){^, n) - Vi{n)(l)li^^){t, n). (5.3) 
At the stretched horizon, the first equation of (5.1) can be written as 

70^0 = kdA^xo - V2^o 



{l,m) ) 



with noise correlation 

{i{l,m){t)i{l' ,m'){t')) = 5ll'5mm'd{iQ{t)iQ{t')) = 5ii'5mm'd^oTHS{t - t') 
= Sll'5mm'THS(t - t'). 



(5.4) 



(5.5) 



5.2 Fluctuation Theorem for Scalar Field in Black Hole Back- 
ground 

From the set of the Langevin equations (5.1), we can construct the corresponding Fokker- 
Planck equation of P(xo, a^i, • • • , xn, vi, - ■ ■ , vn, t) with 2A'" + 1 set of variables; 

TV 



ap = d. 



l^P + Ino P 

7o oxq 7o 



(5.6) 



33 



where we defined U{x) ^ \ Ett [(A-^,)' + Vi {ri)xf\. In this expression, we introduced 
a redundant variable xn-^i for convenience, but eventually set xn+i — 0. This equation 
has a solution describing an equilibrium distribution. 



(5.7) 



A general solution to the Fokker- Planck equation can be formally represented in the path 
integral form as. 



P(X0, Xi, t\x'q, x'i,t^O)(X I Y[ VXkC ^^oTh 



N 



N 



k=0 



(5.8) 



t i=l 



The probability that a trajectory — {{x'q, x'l, • • • , x'^^) {xo, xi, ■ • ■ , xn)} is real- 
ized is given by 



t i=l 



N 



dxi ' 



(5.9) 



On the other hand, the probability that the reversed trajectory F* = {(a;o, • • • , xjv) 
(xq, x'y^ - ■ ■ ^x'^y\ is realized is given by 



N 



p[f;|x] oc 6-^^^^"*^"^°*°-'^"^-"°^^' J] J]5(x, + 



dxi ' 



(5.10) 



t i=l 



Hence the ratio of these two probabilities becomes 



exp 



exp 



/ dtxok{xi - Xo) 



H JV^ 



dt4>\l,rn){^e)drA\l,rn){^e) 



(5.11) 



Here we have written the index for the angular momentum (/, m) explicitly. Summing 
over all the contributions from various partial waves {l,m), the ratio can be written as an 
integral over the stretched horizon; 



n§M^=exp 



{l,m) 



H JVr 



exp 



JTr 



(5.12) 
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Here we have used the definition of the energy-momentum tensor = dt(j)^ (f)'^ — 
dt4>^dr,4>^ ■ Logarithm of the ratio is proportional to the energy fiux into the black hole 
AM[r^] = Jp^ dtdn rlTl[{t, r^, Vi). Hence, by using the first law of black hole thermody- 
namics Tff AjS'Bjjfr,-] = AM[r^], we can interpret this entropy production as an amount 
of difference of the black hole entropy during i = ~ r, 

^£|^ = exp[A5^^M. (5.13) 

In a more general setting, we can introduce an externally controlled parameter such 
as a variable mass term m{t) in the potential U{x; Af ). Even in the presence of such an 
external parameter, the ratio can be shown to be given by the difference of the entropy. 



p^[r*|x] 



exp[A5B^^[r,]]. (5.14) 



In a case with time-dependent external parameters, the forward and the reversed protocols 
are generally different and we need to put F and R to distinguish them. 

In order to apply the fiuctuation theorem, we further multiply the above probabilities 
P^[r,-|x'] (or P^[r*|x]) by probabilities for the initial distributions. As we discussed 
above we can assume that the system is in an equilibrium distribution at the external 
parameter (or Af) with the Hawking temperature P^^{x']\q) (or P^{x;\^)). Hence 



P^[T%\x]P^'i{x] Af) 

- exp [ASBH[Tr] - P {H[x'; A^] - H[x; Af ]) + /3 (F(Ao^) - F(Af ))] 

= exp [{ASbh + ASmWA] ■ (5.15) 

Here, we defined the entropy difference of the matter by A.Sm = —l3{H[x'; Xq] —H[x; Af ])-|- 
/3{F{Xq) — F(Af )), where H[x'; A^] is the total energy of the system at t — with an 
external parameter A^ and F{\q) is the free energy defined by Z{\q) = e~^^^^o). 

The fiuctuation theorem is a direct consequence of the above key relation (5.15). As 
we saw in sec.3, it is straightforward to prove that 

P^{ASbh + A^m) ^ ASbh+^Sm 

Here p^{ASbh+ASm) is the probability to observe a value of the total entropy production 
^Sbh + A5'm with the forwardly controlled external parameter. The denominator is 
similarly defined as the probability to observe a negative value of the entropy production 



35 



in the reversed protocol. Since the right hand side is usually much bigger than 1, the 
numerator is generally much bigger than the denominator. 
By integrating it, we have the Jarzynski equality; 



We observe that there must exist a path with {ASbh + ^Sm) < 0, i.e. an entropy 
decreasing path, otherwise the Jarzynski equality cannot be satisfied. As we saw in 
section 3, the generalized second law [31] 



is derived using the Jensen inequality (e^) > e^^K The above theorems (5.16) and (5.17) 
are also applicable to dynamical processes which are generally in non-equilibrium distri- 
butions, if the Fokker-Planck equation we have used is valid. As we noticed, the validity 
holds when the time scale of the dynamics is longer than the time scale of the inverse 
Hawking temperature h/{2'KTu)- The condition is not always satisfied, and in such situ- 
ations, we need to take effects of time-correlations of emissions. 

5.3 Memory Effect and Quantum Corrections 

In the previous sections, we have approximated the dynamics of the scalar fields by 
Langevin and Fokker-Planck equations. The approximation is valid when the noise cor- 
relation (4.84) can be replaced by the white noise and also the evolution of the scalar 
field (ff is dominated by the classical path described by the Langevin equation. The first 
condition is violated for a shorter time scale than h/27rTH- The second condition is re- 
lated to a justification of the Markovian process we have used. If we take — )■ limit 
while keeping Th = hK/27i fixed, both conditions are satisfied. If these conditions are 
violated, we need to treat the system quantum mechanically without using the classical 
stochastic equations. More detailed analysis of such quantum corrections will be reported 
in a separate paper, and a brief sketch is given here. 

We start from the action (4.78) at the stretched horizon. Before integrating out the 
variable 0", this gives an amplitude of the stretched horizon variables 0^ and 0^. But in 
terms of the variable (f)^' , the path integral represents the evolution of a density matrix 
a la Schwinger-Keldysh, and the path integral (4.79) should be regarded as giving a 
probability, not an amplitude for the configuration 0*". Based on this interpretation, we 
wrote it as P. 




(5.17) 



{{ASbh + ASm))>0. 



(5.18) 
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The classical limit with Th fixed corresponds to replacing K{t,t') by Tf/(5(t — t'). In 
this limit, the probability for a trajectory rT-[0^] with an initial value (pl^^ to be realized 
is given by 



p[r. 



rlV^hii] 



exp 



4T, 



dtdVtrllidt - dr*)(t>''{t)X 



t,r>re,{l,m) 

(5.19) 



The ratio of the forward and the backward probabilities is now given by 



p[r. 



p[r*|0: 



■fin] 



exp 



1 



dtdflrl^'{t)dr*(l)''{t) 



(5.20) 



and reproduces the previous result (5.12). The exponent is proportional to the energy 
flowing into the black hole across the horizon, and interpreted as an entropy increase of 
the black hole. 

More generally, if we do not replace K{t, t') by the white noise, the ratio becomes 



rid] 



p[r* 



exp 



^finJ 



dtdt'dnrl^''{t)K-\t,t')dr*(p^{t') 



(5.21) 



which is nonlocal in time. By expanding the kernel in terms of derivatives of the delta 
functions, the exponent receives corrections to the energy flow. These corrections can 
be interpreted as flows of higher-spin currents (operators containing higher derivatives of 
flelds) into the black hole. These terms vanish after taking a long-time average, but remain 
for a short time scale. Applying the fluctuation theorems with the nonlocal modiflcation 
of the kernel, the entropy increase of the black hole receives higher derivative corrections. 
A geometric interpretation of these corrections is interesting. 

Another important quantum correction is the violation of the Markovian assumption. 
If the path integral is not dominated by classical paths, we need to sum over all possible 
sequences of conflgurations at the level of amplitudes, instead of considering probabilities 
at the classical level. We also need to generalize the fluctuation theorem themselves at 
the fully quantum level. We hope to come back to these issues in near future. 



5.4 Steady State Fluctuation Theorem in Black Holes 

So far, we have applied the fluctuation theorem to a scalar field in an equilibrium distri- 
bution and disturbance around it. In realizing such a situation, we have put the black 
hole in a box with an adiabatic (insulating) wall. Instead we can consider a steady state 
with a constant (but very small) energy fiow from a black hole to outside. This can be 
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realized by putting a black hole in a box in contact with a thermal bath with a slightly 
lower (or higher) temperature than the Hawking temperature. We set the temperature 
at the wall of the box as Tu,{< Th). Since the energy flow is assumed to be very small, 
we neglect backreactions of the energy transfer to the black hole itself. The fluctuation 
theorem for the steady state is reviewed in the appendix C. 

First we start from the following discretized form of the equations of motion; 



7oi;o = -Hxo - xi) - V2^o , {Ut)Ut')) = 7oT//5(t - t') 
xi = -k{xi - X2) + k{xo - xi) - Vi{ri)xi 

Xn = -JXn - k{xN - xn+i) + k{xN-i - xn) - Vi{rN)xN - 

, mm) = iTjit - 1') 

xn+1 = 0. 



(5.22) 



7o k 



Xo 



k k 



Xl 



xn-1 Xn 



7- 



Figure 8: A discretized model of a scalar field in a box. The wall at r = is in contact 
with a thermal bath with temperature slightly lower than the Hawking temperature. 
Then there is an energy flow from the black hole to the outer thermal bath. 



The only difference from the previous model is the equation for the variable a; at at 
the wall. The variable in this model interacts with the thermal bath at temperature 
Tyj. A redundant variable xn+i is introduced for simplifying the above equation. The 
corresponding Fokker-Planck equation is given by 



xo 



7o 0x0 7o 



iV-l 



IVn 



dU 



dxN 

The solution can be written in terms of the following path integral 

Jxi^i J ,1^1 



(5.23) 



X e 



7^ Svr dt^XN+lXN+dj^^uf 



\xn+i=0- 



(5.24) 
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The ratio of the probabihties for a single partial wave with m) is given by 



P[T*\xf] 



exp 



dtxod^gU - 



H JVr 



T 

■J- t 



dtxNixN + d^j^U) 



W JVr 



(5.25) 



|a;jv+i=0 



In addition to the energy flow at the horizon, there is another contribution from the wall. 
The potential U {x) is written as a sum of three terms; 



U{x) = Ui{Xo,Xi, ■ ■ ■ ,Xn-i) + Ui2{xn-1,Xn) + U2{xn,xn+i) 



(5.26) 



where 



Ui{xo,Xi, ■■■ ,Xn- 



1 



i=l 



d 



+ Vi(ri)x'' 



jj , . _l f Xn - Xn-1 

Ui2{Xn-1, Xn) — - 



U2{xpf, Xn+i) 



d 



Xn+1 — Xn 



d 



+ Vi{rN)x% 



(5.27) 



We turn on the potential U12 at the wall during a time interval between t = and 
t — T. This can be realized by introducing the external parameter controlling the potential 
U12 such as 



Ui2{xN-\,XN-At) = ^ (V ~ 1^ ~ il) Ui2{xn-1,Xn)- 



(5.28) 



Then the variables {xq^xi,--- ,xn-i) are decoupled from xn when t < and t > t. 
Since the external thermal bath is decoupled for a long time during t < 0, the state can 
be considered in the equilibrium at t = 0. The ratio of the probabilities of the initial 
distributions is, hence, given by 

1 /l"^"' 



pcq_(rf.. ^ 



P^{Xi,n) 



exp 



H 



1 



2 (*iV,ini - ^Nm) + ^2(2;ini) " U2{,Xf,^) 

The second terms are canceled by the following terms in eq.(5.25) 

"1 



(5.29) 



dtXNixN + dx^U2{xN)) 



-mxlf + U2{xn) 



fin 



(5.30) 



The remaining terms in (5.29) is, of course, independent of the duration r, and can be 
neglected compared to other terms in (5.25) that are proportional to r. 
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As a result, if take the leading contributions in the large r hmit and neglect C(r°) 
terms in the exponent, we obtain 



-f l-^ini] ''(•^ini) 

P[r*|xfin]P^q(a;fi„) 



exp 



exp 



1 



dtxod^^Ui - 



1 



T 



dtX]sfdxf^Ui2 



dtxokdA Xq 



T 



dtxokdA Xn 



w JT 



(5.31) 



Because of the energy conservation for a steady state, we have the relation J dtxNA^XN = 
— J dtxoA~XQ. In the contimuum limit N ^ oo with the scalings explained before, the 
logarithm of the above ratio becomes 



dtxokdA Xq 



H Jr-r 



T 



dtxNkdA Xn 



W JT-r 



{/3w - Ph) / dtdt(t/n^){t, r,)dr,(pa^^){t, r^) = A/3r J(,,^)[r, 



(5.32) 



We have defined = — Ph, which is positive from our assumption Th > Ty,. By 
summing all the contributions from the partial waves with {l,m), we have 



J[rr] ^ 



dtdQriTl{t,r,,n), (5.33) 



where we have used the definition of the energy momentum tensor T[ — dtcj)''' d""' (jf — 
dt(t>^ dr^(tf . J[rT-] is a current flowing at the horizon out of the black hole. Prom the setting 
Th > Tu,, the averaged current is positive, but it can take both positive or negative values 
because of fluctuations of absorption and emission by the Hawking radiation. 

We now have established the steady state fluctuation theorem in the black hole back- 
ground as 



lim — 

T—^OO T 



= A/3Joo. 



(5.34) 



For the deflnitions of p, see eq. C.6. The theorem can be restated in terms of a gener- 
ating function Z{ar, Af3), and leads to various relations between the response coefficients 
L(i), L(2), . . . deflned by (J^) = L^^^Ap + L(2)/2(A^)2 + • • • and correlator of currents 
{J{t)J{t')) . For more details, see the appendix C. 

In our case, these relations lead to the following relations; 

dt I dn rt{T[{t, r,)T[{0, r,))|A/3=o (5.35) 



L(2) = ^lim i- fdtdt'J dil rtdAp{T[{t,r,)T[{t\r,))\Ap=o (5.36) 



40 



The first relation is the Green- Kubo relation for the energy current fiowing at the horizon 
r = Tg. The second one is a non-linear generalization, and the evaluation of the right 
hand side needs the derivative of the correlation function with respect to the temperature 
difference. This means that the information of the equilibrium distribution at = is 
not sufficient to obtain the non- linear response function of the current. 

6 Summary 

In this paper, we derived a stochastic equation with a dissipative term and a noise for a 
scalar field in a black hole background. The dissipation comes from the ingoing boundary 
condition at the horizon while the noise comes from the Hawking radiation. The stochastic 
equation can be derived by considering a stretched horizon and integrating variables 
between the horizon and the stretched horizon. The stochastic equation describes the 
dynamics of the scalar field in the limit H ^ with the Hawking temperature Th — Hk /27r 
kept finite. We then applied the non-equilibrium fiuctuation theorems, developed in the 
statistical physics, to the above stochastic equation in the black hole background. We 
consider two cases. One is a scalar field confined in a box with an insulating wall. The 
system is relaxed to an equilibrium state at the Hawking temperature. The fiuctuation 
theorem leads to the second law of the black hole thermodynamics after taking a thermal 
average. The other case is a scalar field in a box in contact with a heat bath with a 
different temperature from Th- Then there is an energy flow between the horizon and 
the outer boundary. The fiuctuation theorem leads to the Green-Kubo relation and its 
nonlinear generalizations. 

We have used an approximation of replacing a nonlocal (colored) noise correlation by a 
white noise. We furthermore approximated the dynamical evolution of the scalar field by 
a classical Markovian process. These approximations are valid in the classical limit h ^ 
with the Hawking temperature Th flxed. In this sense, quantum effect is partially taken 
into account through the Hawking radiation. The results such as the ordinary second 
law of the black hole thermodynamics or the Green-Kubo relation are derived only in 
such approximations. As mentioned in Sec 5.3, the nonlocal noise correlation leads to a 
deviation of the black hole entropy appearing in the second law of thermodynamics. We 
hope to come back to these issues in near future. 

Finally it is interesting to generalize our results for a scalar field to a vector or a 
gravitational field and obtain quantum corrections to the membrane action [4, 32]. The 
absorption of energy across the stretched horizon is known to give dissipative equations 
such as the Ohm's law or Navier-Stokes equation on the membrane. If the Hawking 



41 



radiation is included, these equations must receive quantum corrections as noise terms. 
Then the Hawking radiation may be interpreted as the anomaly inflow of the gravitational 
and gauge anomalies. In this sense, the quantum membrane action will be analogous to 
the edge state of quantum Hall effect. It is furthermore interesting if we can relate such 
a quantum membrane interpretation to the black hole entropy. 
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A Path integral form of the Fokker-Planck equation 



In this appendix we derive the path integral form (2.11) of the solution to the Fokker- 
Planck equation; 



dtP{x, V, t\xQ, Vq, 0) = LptP{x, v,t\xQ, Vq, 0) 
= i-vP) + 
For a small time- interval, it can be written as 



— + — ^ 1 P 

m m ox 



+ a,^(gp). (A.l) 



P{x, V, At\xo, vo, 0) = e^^^^''S{x - xo)5{v - vq) 



72^ 



1 + At 



^ikx{x-xo)+ikv{v-vo) 



xexp[-^(mH^ + 7.o + ^)' 



(A.2) 
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Then by using the Chapman- Kolmogorov equation P(X3\Xi) = J (iX2P(X3|X2)P(X2|Xi) 
which is equivalent to an insertion of the complete set and integrating over v, we obtain 
the path integral form as follows; 



P{x,t\xo 



J X 



x{t)=x 



Vxexp 



x{0)=xo 



-4^ / dt'{mx + ^x + %y 



(A.3) 



If we use the Langevin equation (2.1), the path integral is equivalent to the noise average 
with the weight function in eq. (2.3). 



B Noise correlation and Hawking radiation 

The noise correlation induced in the effective equation of motion for the boundary field 
at the stretched horizon r — rn + € can be interpreted as the Hawking radiation. Here 
we first review the method to determine the energy-momentum tensor in the black hole 
background by using the trace anomaly of the energy-momentum tensor and the regularity 
condition at the horizon [33], and then generalize the method to determine higher spin 
currents [34, 35, 36, 37]. 

In two dimensions, the trace of the energy-momentum tensor of a single scalar field 
(i.e. the central charge is c = 1) has an anomaly term proportion to the scalar curvature 
R 



247r 



R. 



(B.l) 



Writing the metric in the conformal gauge ds^ — e'^^'"''"'{—dudv), the equation becomes 
^dudy(f. By combining with the conservation of the energy-momentum tensor 



T 

-*- II.- 
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V^T^ = 0, derivatives of the EM tensor dyTuu{u,v) and du%y{u,v) can be written as 
follows; 



d T 



d T 



247r 
1 

24^ 



Prom these equations, we can define a (anti-) holomorphic quantity 



(B.2) 
(B.3) 
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(B.4) 
(B.5) 
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They are often called (anti-) holomorphic energy-momentum tensors, but their trans- 
formation laws are anomalous and not tensors in the exact sense. Actually, under a 
coordinate transformation from {u,v) to {U = U{u), V = V{v)), they transform as 



tuu{U) 



[kUY 



(B.6) 



where {U, u} is the Schwarzian derivative. 



(B.7) 



In particularly, for the transformation from the Schwarzschild coordinates to the Kruskal 
ones, namely from {u,v) to {U,V) = {—K~^e~'^'^,K^^e'^'"), the Schwarzian derivative be- 
comes {U,u} = 

Now, we impose the regularity condition at the horizon. The energy momentum tensor 
Tjju must behave regularly near the future horizon = in the regular coordinates, and 
so is tuu{U) since they are related regularly as (B.4). The regularity condition, hence, 
imposes that tuu must behave as 

tuu{u ^oo) = — . (B.8) 

4o7r 

If we neglect the effect of scatterings of the outgoing fluxes (namely in the absence of the 
gray body factor), we can extrapolate the above flux at the horizon to the outgoing flux 
at r ^ oo as 

,2 



T„„(r oo) 



K 



H- 



487r 12 

It is interpreted as the flux from the black body with the Hawking temperature T^, 



(B.9) 
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The transformation property of the holomorphic energy-momentum tensor can be also 
derived by considering the following point-splitting regularization. 



: tun{u) : = lim [ducl>{u + '^)d^ct>{u - f ) - {duct>{u + '^)d^ct>{u - f ))] 

1 



lim 



a„0(M + |)a„0(«-|) + 



47r(52 



(B.ll) 



where we have used the explicit form of the free boson propagator {(j){u)(j){u')) — — \a{u — 
u')/A7r. Prom this deflnition, we can relate it to the energy momentum tensor regularlized 
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in the Kruskal {U) coordinate; 



lim 

5^0 



lim 

5^0 



d^U{u + l)d^U{u - l)du<P{U{u + l))duHU{u - f )) + 4^ 



duU{u + l)duU{u - I) tuu{U) - 



47r ([/(m + I) - - |))2 / 47r52 



{duUf : tuu{U) -.K-T-riU^u} 



247r 



(B.12) 



Namely, the Schwarzian derivative is nothing but the difference of the normal orderings 
in different coordinates. 

The energy fiux (which corresponds to the first moment of the thermal spectrum 
(B.IO)) can be generalized to a fiux of a higher spin current with a higher moment, and 
its generating function can be defined as a correlation function of the scalar field; 



J(u, u + a) = J2-[- 5u</'(m)5"+V(«) : 
=: du(t>{u)du(l){u ^ a) : . 
The normal ordering : : is defined similarly to tuu{u) by 

: du^{u)du^{u) : = lim [du^{u)du^{u') - {d^^{u)du4>{u'))] 



(B.13) 



(B.14) 



Then we can show that J(m, u + a) transforms under the coordinate transformation from 
uioU {u) as 



J{u, u + a)^ duU{u)duU{u + a)J{U{u), U{u + a)) + 



K 



4 sinh^ ^ ' 



+ 



(B.15) 



Similarly to the energy flux discussed before, the regularity condition at the future horizon 
fixes the value of J{u, w + a) at C/ = as 



J{u, u^-a) \r=rH du(f){u)du(j){u + a) := 



47r 



K, 



4 sinh 



2 «a „2 



1 



(B.16) 



This can be interpreted as a correlation function of du4>{u) and du(f){u + a) on the Kruskal 
vacuum. 

In Section 4.3.3, we have shown that the scalar field obeys a stochastic equation of 
motion 



du(l>{t-r*)\r=r„+e^-V2^{t). 



(B.17) 
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at the stretched horizon. Since the equation is independent of the value of e, we can safely 
take e — >■ limit. Then the value of the generating function J{u, u + a) for the higher 
spin fluxes discussed above is equivalent to the noise correlation 2{^{t)^{t + a)) of the 
Langevin equation at the horizon. The functional forms are equal, though the coefficients 
are different by a factor 4. Reason for the discrepancy is now under study. 



C The Steady State Fluctuation Theorem 

In this appendix, we consider the fluctuation theorem for a steady state and derive the 
Green-Kubo formula. 

Assume that we have two variable Xi,X2, and each of them is in contact with a different 
thermal bath with temperature Ti and T2. We further assume that the dynamics is 
governed by the set of Langevin equations such as 

dV 

miVi + — + ^,v,^Ci , {Cimi{t'))^'2^iTiS{t-t') 
dV 

m2V2 + ^ + 72V2 = ^2 , (6(^)6(0) = 272^25(^-0- (C.l) 
0x2 

Here, 2:2; Af ) is an interaction potential between the two variables. The corre- 

sponding Fokker-Planck equation of the system can be obtained straightforwardly. The 
trajectory T^- is also generalized as = {^i^) = i^iit): ^2{t))\x{0) = (a;i(0), 0:2(0)) = 
(Xj^j, xfnj)}. Then the solution to the Fokker-Planck equation gives probabilities of the 
forward and the reversed protocols, and the ratio is given by 



P [rT-|Xini] 



1 /■ , . / .. dV{x;\f)\ I r , . { .. dV{x;\f] 
dtxi rriiXi H — — / atx2 777.2X2 -l- 



Ti Jt, \ dxi J T2 Jr^ V 9x2 

(C.2) 



We have assumed that the two variables are decoupled before t — and after t — t; 
the interaction potential V vanishes at t < and t > r. The initial distribution of 
the total system is given by a product of the equilibrium distributions of each system 
P^'^ixini) — -P^'^(2^ini)-P^*^(3^hii)- The forward protocol is expressed as 

V{x- Af ) = Vi{xi) + V2{x2) + f{Xf)Vu{xi - X2) (C.3) 

where 

/(Af) = ^(y-|Af-^|), Af = ^. (C.4) 
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T_ means r — e for < e <C r. Function f{t) satisfies f{t = 0) = f{t = r) ~ and 
/(O < \t\ < t) = 1, so that the interaction switches on at i = and off at t = r. This 
protocol has the reversal symmetry /(Af ) = /(Af_J. 

In considering the large interval limit r — > oo, the energy transfer such as J dtxidx^Vi2{x\ — 
X2) (or J dii;29j;2T42(3:^i— 2^2)) grows linearly in T. On the other hand AE'i = J dtxi{miXi-\- 
dx^Vi{xi)) = {\rnixl+Vi{xi))t=T-{\rnixl+Vi{xi))t=o or AE^a = / dtx2{m2X2+dx^V2{x2)) 
is at most C(t°). If each system becomes stationary after taking r — )■ 00, the change in 
the energy of each system vanishes. Hence we can drop the term P^{xird) / P^ixf^ry) and 
a contribution of AE'j in P[rT-|xini]/i-*[r*|xfin] when we evaluate the quantity 

Hm-mf^l^^N^^V (C.5) 
™t \P[V%\x^^]P^^{x^r) ) ^ ' 

In addition, we have dtxidiVi2 ~ — /r dtx2d2Vi2 + ^(t^). Therefore we can write 
the ratio of the probabilities only in terms of the energy current defined by J\^t\ = 
^ Jp dtx idiVi2. Writing the temperature difference as A/3 = /32 — we obtain the 
following relation; 

p(J„ A/3) = j Pa;P[r,|Xini]P^^(Xini)5(J-r - J[r,]) 

~ j Pa;P[r;|a;fijP^^(a:fin)e-^^-^"P-]5(J,- J[r,]) 

= e-'^^'^ j VxP[Ti\x^^]p''\x^^)8{Jr + J[r;]) 

^e'^^'^p{-Jr.m. (C.6) 

The steady state fluctuation theorem can be written as 

p(J,,A/3) 



lim — In 

r— >-oo T 



From this relation, we can derive the Green-Kubo relation and its non-linear general- 
izations. By using the generating function 

Z{a,, AP) = In (^y" dJre'-^'^'^^piJr, A^)^ , (C.8) 

the steady state fluctuation theorem (C.6) can be recast into 

Z(ar + iAp,A/3) ^ Z(-ar,A/3). (C.9) 
Taking a derivative of both sides with respect to A/3 and setting A/3 = 0, we have 

9a/3 [Z{ar, 0) - Z{-ar, 0)] = -ida,Z{ar, 0). (C.IO) 
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The generating function can be expanded in terms of the correlators of Jr as 

A^) = ^^^G„(A^). (C.ll) 



n=l 



Gn{/3) gives a connected Green function of the averaged current 



Jr = - f dtJ{t). (C.12) 



T _ ^ 

Now the equation (C.IO) is rewritten in the following form; 

[1 - (-1)"] aA/jGn(0) = TGr,+M- (CIS) 

We further expand the one-point function of J,-, which gives an expectation value of the 
current, with respect to the inverse temperature difference as 

°° 7-(m) 

Gi(A/3)^ J]— -(A/3r. (C.14) 
m! 

For n = 0, wc have a trivial identity Gi(0) = L'^'^^ = 0. For n = 1, the Green-Kubo 
relation is derived; 

L^'^ = ^ f dtdt'{J{t)J{t%Ap=o 

dt{J{t)JmiA0=o. (C.15) 

When A/3 = 0, the system is described by the equilibrium distribution function P^'^{x) = 
g-/3£;tot(a;)^2', /3 = /^i = and an expectation value of a function F[x{t)) is given by 

{F{x{t)))\Ai3=o = J VxP^\x{t))F{x{t)). In the larger limit, the correlator (J(t)J(t')) I A/3=o 
depends only on [t — t'). 

We can also obtain the expression of L^^^ , L^^^ , • • • by taking further derivatives of the 
equation (C.9) with respect to A^fl. For instance, we can derive 

d% [Z{a^, 0) - Z{-a^, 0)] = -id^^dAp [Z{a^, 0) + 0)] 
^ (1 - (-1)'^) ai^i^„(0) = T (1 + (-1)"+^) dApK^+m. (C.16) 

For n = 1, we have 

L(2) = hm ^ /" dtdt'dAf){J{t)J{t'))\A^=o. (C.17) 

These nonhnear generalizations can be systematically obtained by using the steady state 
fluctuation theorems. We apply these expansion method to a system of a black hole and 
matter to obtain the Green-Kubo relation for a thermal current in the body of the paper. 
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